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We study order-parameter fluctuations (OPF) in disordered systems by considering the behavior 
of some recently introduced paramaters G, G c which have proven very useful to locate phase tran- 
sitions. We prove that both parameters G (for disconnected overlap disorder averages) and G c (for 
connected disorder averages) take the respective universal values 1/3 and 13/31 in the T —* limit 
for any finite volume provided the ground state is unique and there is no gap in the ground state 
local-field distributions, conditions which are met in generic spin-glass models with continuous cou- 
plings and no gap at zero coupling. This makes G, G c ideal parameters to locate phase transitions 
in disordered systems much alike the Binder cumulant is for ordered systems. We check our results 
by exactly computing OPF in a simple example of uncoupled spins in the presence of random fields 
and the one-dimensional Ising spin glass. At finite temperatures, we discuss in which conditions the 
value 1/3 for G may be recovered by conjecturing different scenarios depending on whether OPF 
are finite or vanish in the infinite-volume limit. In particular, we discuss replica equivalence and its 
natural consequence limv^oo G(V,T) = 1/3 when OPF are finite. As an example of a model where 
OPF vanish and replica equivalence does not give information about G we study the Sherrington- 
Kirkpatrick spherical spin-glass model by doing numerical simulations for small sizes. Again we find 
results compatible with G = 1/3 in the spin-glass phase. 



I. INTRODUCTION 

It is well known that mean-field spin glasses are characterized by strong (non vanishing in the thermodynamic limit) 
sample to sample fluctuations of the order parameterEJ. Despite the fact that extensive thermodynamic quantities 
(such as free energy and all its finite order derivatives) are self-averaging in the thermodynamic limit (i.e. their 
intensive part does not depend on the realization of the quenched randomness) the same result cannot be extended to 
order parameter fluctuations. It is widely believed that absence of self-averageness of the order parameter is strongly 
related to replica symmetry breaking, i.e. the existence of several ergodic components not related by any symmetry 
of the Hamiltonian. n 

Recently, Guerra suggested^ that sample to sample fluctuations of the order parameter (hereafter referred to as 
OPF) verify some sum rules which are generally valid in any disordered system. This claim assumes that the system is 
stochastically stable in the presence of a mean-field perturbation, a property which may strongly depend on the nature 
of the equilibrium state. A system is stochastically stable if its properties (static or dynamic) smoothly change in 
the presence of a small random perturbation. These sum rules have been recently used to define a new dimensionless 
parameter (hereafter called G) which takes into account sample to sample fluctuations!! This parameter has been 
shown to provide an alternative and powerful way to locate phase transition points in disordered systems. The 
advantage of G respect to more canonical ones (such as the Binder cumulant ratio used in ordered systems) relies 
on the fact that it works very well also in the absence of time-reversal symmetry in the Hamiltonian or.-ather more 
complex disordered systems. In particular, the method has been recently applied for Ising spin glassesBd, Migdal- 
Kadanoff spin glassesQ, Potts glassesQ, Heisenberg spin glasses, which display a chiral phase transitions as well as 
some protein folding models! . 

The purpose of this paper is to show, by using general arguments, analytic computations for simple models and 
numerical simulations, that indeed this new parameter is the appropiate tool to investigate phase transitions in 
disordered systems much like the Binder cumulant is for ordered systems. We conjecture and prove that this parameter 
G takes the universal value 1/3 at zero temperature for any disordered system (finite or infinite) with the only condition 
of uniqueness of the ground state and absence of a zero-temperature gap in the local field distribution. This condition 
is satisfied by all spin-glass models with continuous distribution of couplings and no gap at zero coupling. At finite 
temperature G certainly depends on the system size. We claim that due to the property of replica equivalence, for 
models in which OPF are finite, G converges in the infinite-volume limit to zero if the system is in a paramagnetic 
phase and to the same zero-temperature value 1/3 if the system is in the spin-glass phase. When OPF vanish this 
does not necessarily hold and we discuss in what conditions the universal value 1/3 may be recovered. 
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The paper is organized as follows. Section II is a reminder of the definition of the G parameter as well as some other 
useful ones. Section III presents a detailed computation on a simple disordered model which serves as an illustrative 
example of the main results. Section IV proves the zero-temperature conjecture under some general conditions for any 
disordered system. Section V presents detailed calculations on the one-dimensional Ising spin-glass model using the 
transfer matrix approach. Section VI addresses the validity of the conjecture at finite temperature by studying the 
Sherrington-Kirkpatrick spherical spin glass, a model where OPF vanish. Finally we discuss the results and present 
the conclusions. 



II. THE G PARAMETER AND REPLICA EQUIVALENCE 

The definition of the G parameter is based on some exact relations, obtained for the sample to sample fluctuations 
of the order parameter in the Sherrington-Kirkpatrick (SK) modeffl. The SK model is defined by the disordered 
mean-field Hamiltonian, 

Ti-SK = ~ / ] Jjj GjVj > (1) 

i<j 

where the are quenched Gaussian variables with zero average and variance 1/N where N is the number of sites. 
The SK model presents a second order phase transition at T e = 1 below which replica symmetry breaks down and 



ergodicity is broken. The spin-glass phase is described by an order parameter function Pj{qi2) where q\2 = $2»=i °\ a \ 
is the replica overlap and the subindex J stands for the realization of the quenched randomness. Pj{q) is a simple 
object in the paramagnetic phase above T c (Pj(q) = 5(q)) but develops strong sample to sample fluctuations below 
inside the spin-glass phase. Fluctuations-in the order parameter were originally derived by Bray, Moore and Younj 
using the property of replica equivalence^. This property states that the sum of all elements contained in a given row 
(or column) in the replica matrix Q a \, is independent of the row (or column). As shown by Parisilij this is a necessary 
condition for the replicated free energy to be proportional to the number of replicas n and have a well defined free 
energy in the limit n — > 0. Fluctuations are then described by the following exact relation in the N — > oo limitH, 



Pj(qi2,q3i) = ^P.j{qi2)5{qi2 - q 3 i) + ^Pj{qi2)Pj{qu) , (2) 

where (.) stands for disorder average and 1,2,3,4 denote replica indices. Therefore, 



Pj{qi2,qu)^P.j{qi2)P.j{qu) , (3) 

so Pj fluctuates with in a non-trivial way. Multiplying both sides of ca-ph by q\ 2 and q^ 4 and integrating over 
all possible values of the overlaps qi 2 , (734 one obtains the following sum ruleEHl, 

J&=lW) + l~(^) 2 ■ (4) 

where (...) stands for thermal average. This relationship has been also rederived by Guenja using general arguments 
based on self-averaging properties of the internal energy as well as its finite derivatives^. Now let us define the 
following ratio, 
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G = ^ ' g n (5) 

<9 4 > - (<Z 2 > 

Note that the numerator in (^) corresponds (except for the absence of a multiplicative constant iV 2 ) to the sample 
fluctuations of the spin-glass susceptibility. For the SK model, because of the sum rule (H) , it is possible to show that 
G takes only two values. G is equal to 1/3 in the replica symmetry broken phase and vanishes above T c , 

G = ie(T c - T) . (6) 

The generality of the replica-equivalence property suggests that (g) will hold in any system (even beyond mean- 
field) if OPF do not vanish in the limit V ■— > 00. But may well happen that OPF vanish. Then both numerator and 
denominator in ([5]) vanish in the V — + 00 limit. In this case replica equivalence is not enough to decide what the value 
of G is. The value of G is then determined by the form of the finite-size corrections to the order parameter (and in 
particular its prefactors), which in principle could not satisfy sum rules such as (^|). Despite this uncertainty, in this 
paper we propose three possible scenarios for the parameter G, 
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• 1. OPF remain finite in the thermodynamic limit. This is the general situation encountered in mean-field models 
with a replica broken phase. So both numerator and denominator in (||) are finite in the infinite volume limit. 
The property of replica equivalence and also stochastic stability indicate that the same should be valid for any 
finite-dimensional disordered system (assuming that for those systems OPF are finite) leading to G = 1/3 in 
the spin-glass phase. 

• 2. OPF vanish in the large volume limit like 1/V. This is the situation typically encountered in the paramagnetic 
phase. The ratio may then be zero or finite depending on the particular case. 

• 3. OPF vanish in the large volume limit slower than 1/V (for instance, like y% with a < 1). This situation is 
typical of disordered systems with a marginally stable replica symmetric phase. Both numerator and denomi- 
nator in ([5]) vanish, the ratio G is finite but may be different from 1/3 at finite temperature. In this case the 
property of replica equivalence cannot be used for the reason discussed before and stochastic stability may not 
hold. Actually the property of stochastic stability may breakdown if the equilibrium phase is drastically changed 
in the presence of a mean-field perturbation. This situation may be found in spin-glass models without OPF such 
as hierarchichal lattices (i.e. spin glasses in the Migdal-Kadanoff approximation), the Sherrington-Kirkpatrick 
spherical spin glass (see section VI) or finite-dimensional models described by a unique low-temperature state 
such as the droplet model. 

Despite the main hypothesis of stochastic stability remains still to be proven all previous three cases seem quite 
reasonable and we do not know of non-trivial counterexamples. Note that there is not any direct relationship between 
OPF and the value of G in the low-temperature phase. Actually, the previous possibilities 1 and 3 may yield the samie 
value of G although the physical description of the low-temperature phase is much different. As has been observed incl 
the non vanishing of G should not be taken as direct evidence for non-vanishing OPF or replica symmetry breaking. 
In order to better evidenciate whether OPF survive in the infinite- volume limit, it is necessary to consider another 
dimensionless parameter, which has not the ambiguity of the ratio of two vanishing quantities. For instance one may 
define the A parameter^, 

A= w-F> 2 , ( 7) 

which is nothing else than the numerator of (^|) appropriately normalized. We will see later that the nice properties 
of G are not present in the parameter A and the former is much more convenient to locate phase transitions. Generally, 
one expects A to be a non trivial function of both volume and temperature vanishing (in the V — > 00 limit) only when 
OPF vanish (for instance, in a paramagnetic phase). If OPF are finite A may take a finite value but an identity such 
as @ for A does not hold. 

In this paper we will show examples for all three behaviors, by explicit analytic computations and some numerical 
calculations. Furthermore, we will show that for models with a unique ground state and without gap in the ground- 
state local field distribution, 

Hm G(V,T) = ± , (8) 

so the G parameter is 1/3 at T = for any finite volume V. This is not anymore true at finite temperature where 
the parameter G may take the value 1/3 only in the infinite volume limit. 

Before finishing this section let us remind that in referencestffl another quantities similar to (^) and (^) have been 
proposed for systems without time-reversal symmetry. These arc defined by considering the connected overlaps, 



q = ((g-(g)) 2 ) 2 -((g-(g)) 2 ) (9) 

(( q -( q m-(( q -( q )r) 2 



Ac = <( g -fa )) 2 ) a -<( g - w)_ . (10) 

((q-(q)) 2 ) 



We will see that a result like (H) applies also to the parameter G c and our result reads: 



HmG (V,T) = H . (11) 
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For the SK model the quantity G c is defined by restricting thermal averages to the q > part of the P(q). G c does 
not satisfy the identity (|6j) so this is not the best quantity to look at in numerical simulations despite the fact that 
both G c and G (and also A c and A) may take similar values in the vicinity of the critical region. This explains why 
similar results were obtained for both sets of quantities in numerical simulations. 



III. AN INSTRUCTIVE EXAMPLE 



Here we analyze in detail a solvable case which will be useful to illustrate the main contents of the paper and how 
disorder expectation values of the overlaps are computed. Moreover, the analysis of this section will prove to be useful 
for a constructive proof of the zero-temperature results @ and ([ll]) to be presented later on. Consider the following 
Hamiltonian, 

v 

H = - s £ j K i a i , (12) 
i=i 

where the spins may take the values ±1 and the fields hi are uncorrelated and randomly taken from a distribution 
P(h) with finite weight at zero field (i.e. -P(O) finite). In principle P{h) may be any function 

P(h)=P{h) + J2 c k6(h-h k ) (13) 

k 

with P(h) any continuous function, all ^ and P(0) ^ 0. This condition is enough to ensure the non-degeneracy 
of the ground state because there is a single configuration which minimizes the energy a* = sign(/ii). Note that if a 
finite fraction of the fields hi were zero then the ground state would be degenerate. With this very general condition 
we may exactly compute the parameter s G and A introduced in the previous section. Everything reduces to compute 
the three overlap quantities: {q 2 ),(q i ) and {q 2 } 2 . The computations are quite elementary and here we present the 
final results. For the numerator and denominator of eq.(|^) we get, 

Numerator = Tr „ — -R 2 + — -R 2 R ^ — -R , (14) 

yo yo yo 

Denominator = - + ^if - 2<y2V -J - + ^ i? 4 , (15) 



where 



R= I dhP(h)t&nh 2 ((3h) , (16) 



R 2 = / dhP(h)tanh 4 (ph) , (17) 



and P(h) is the generic distribution (|13|). The expressions for the parameters G and A may be further simplified 
yielding, 

c _ (W-R 2 )(W+(2V-3)7f) 
(l-i? 2 )(l + (2V-3)i? 2 ) 

and 

A = 2(V -1)(W- gKg + (2V - 3)if ) 
V{l + {V ~l)R 2 ) 2 

Note that in the limit V — > oo both numerator and denominators in ( pT[ ) and ( [l5| ) vanish. The quantity A also 
vanishes like 1/V but the ratio G stays finite, 
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lim G(V, T) = 



R 2 -R 
1-R 2 



(20) 



the asymptotic low-temperature behavior of R and R 2 is given by, 



The finite volume quantity G(T, V) in (18) satisfies the conjecture (g). A simple integration by parts reveals that 

(21) 



R=l — TD + 0{T 2 ) 
where D is a positive constant given by 

D = 2P(0) 



T 



R 2 = R — —D + 0(T ) 



(22) 



Substituting the asymptotic behavior ( pl| ) in ( [l~8"|) we obtain G(V, T = 0) = 1/3. Note that the same result is obtained 
substituting ( pl| ) in ( pp| ) because in this simple example the two limits T — > and V — > oo may be interchanged. 
This is not generally true: in particular when a phase transition takes place at T = the two limits may not be 
interchanged anymore. 

For the parameters G c and A c introduce in (p|)(^0|) we get: 



G c 



R i - (R 2 ) 2 



2V - 2 - 4(V - 2)R 2 + (2V - 3)(i? 2 ) 2 - 3i? 4 



_ fl 4 - (R 2 ) 2 
c V(l - R 2 ) 2 



(23) 



(24) 
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FIG. 1. Parameter G for V = 2, 4, 8, 16, 50, 100 from below to above 
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FIG. 2. Parameter G c for V = 2, 4, 8, 16, 50 from below to above 

We observe that G c behaves in a different way. It tends to zero for T hnite and is again independent of the 
volume for T = but takes the value 13/31. For G c the two limits (V — ► oo and T — > 0) now do not commute. 
Figures 0, |[ [| [l| show the behavior of G, G c , A, A c as functions of temperature for different values of V for the case 
of a Gaussian fields distribution P(h) = (2II) -1 / 2 cxp(— h 2 /2). 
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FIG. 3. Parameter A for V = 2, 4, 8, 16, 50 from up to down 
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FIG. 4. Parameter A c for V = 2 (continuous), 4 (long-dashed), 16 (short-dashed), 50 (dotted) 



In case of a gap of amplitude A in the field distribution one finds that both A and G vanish exponentially with 
that gap G ~ Texp(— (3 A) and the conjecture does not hold anymore. 

We will now prove that, under some general conditions, the conjectured zero-temperature values for G and G c hold 
for any disordered system. 

IV. A PROOF OF THE CONJECTURE 

To generally prove (||) and ( pd| ) we start by considering a general Hamiltonian Ti({a}) where the {af, i — 1, .., V} 
are Ising variables which may take the values ±ll!3. This Hamiltonian may be written in terms of the local fields, 



H 



2_j hid 



1=1 



where the hi are local fields proportional to 



(25) 



(26) 



which depend on the configuration {a}. Suppose now that the Hamiltonian H may only take continuous values so 
there is no ground state degeneracy (apart from a global symmetry in the Hamiltonian such as time-reversal symmetry; 
this case will be discussed below). In particular, no local field hi vanishes. Let us denote by {tr*} the (unique) ground 
state configuration. The ground state is stable with respect all possible number of spin flips, so that the value of 
the energy in that configuration Ti({a*}) is an absolute minimum. In particular the ground state is stable respect to 
single spin flips and the local fields evaluated at the ground state satisfy the property, 



signih*) 

where h* are evaluated at a*. To prove the conjecture we need to prove the following statement: 



(27) 



• Statement: Excitations which involve the reversal of a single spin yield the dominant contribution to the low- 
temperature behavior for all the quantities (q k ) 1 for any positive integers k, I and by extension, to the numerator 
and denominator in (p|,[7p|JlO|). 
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This statement somehow allows to map the most probable excitations in (pq) with those of the instructive example 



presented before (12). Nevertheless, we must emphasize two points. The first one is that the ground-state local field 
distribution in the previous example ([l2]) was taken uncorrelated for different sites and also the same distribution was 
taken for each spin i. In general this is not true. Local fields at different sites may be correlated and the distribution 
on a given site may depend on the site. For instance, in models with open boundaries the local field distribution for 
the sites located on the surface is certainly different from the distribution of those in the bulk. The second observation 
is that, in general, the lowest excitations in ( pof ) may involve groups of several spins (and not a single spin flip like 
in the simple case (|l2|)). So in order to prove the conjecture we must show that excitations in ( pq ) which involve 
the reversal of any number of spins larger than one always yield subdominant low-temperature corrections to the 
single-spin excitation case. 

In what follows we present a constructive proof of the previous statement without need to refer to the results of the 
instructive example which had some restrictive assumptions. We start from the general Hamiltonian ( p5| ) and analyze 
the low-temperature behavior of the order parameters (q 2 ), (q 4 ) and (q 2 ) 2 . We will first consider the case of one spin 
excitations and later on the more general one of higher-order excitations. 

• One-spin excitations 

The calculation proceeds as follows. Consider the ground state {a*} of (|2^) as unique and one-spin excitations 
which involves the reversal of a single spin. If we consider the ground state plus this class of V possible excitations 
we can compute the correlation function (o~iO~j) (i ^ j), obtaining the result: 

/ exp(- 2(3h* a*) + exp(-2/?/i*cr*) \ 
(a t a 3 ) = a*a*{\ - 2 ^ v 3 J ) (28) 

Defining x; t = cxp(—2(3h*a*), we get in the f3 — > oo limit, 

(Wi) = or***(l - 2( Xl + x 3 )) (i + j), (29) 



where we have aproximated by 1 the term in the denominator of the ratio in (|28|). Such an approximation is 
allowed provided one performs the limit (3 — > oo before the infinite volume limit. Note that, in that denominator, 
each exponential contributes to the sum at most with a term proportional to the temperature (see below). 
Because there are V terms of that type, at most that term is of order VT. Hence, in the limit TV <C 1, that 
denominator equals 1. The result ( p9| ) is the only quantity we need in order to evaluate (q 2 ) and (q 2 ) 2 - In terms 
of the variable = (cjCj) 2 , these are given by, 

<S 2 > = ^ + ^f£^ . (30) 

<9 2 > 2 = ^ + 73E^ + |iE^ + ^ E T ^ + ^z E T ^ > ( 31 ) 

where the indexes in the sums run from 1 to V and correspond to different sites. To average ([50|) , (||l]) over 
the disorder we need to compute disorder averages of terms of the type x™x™ where i,j denote sites and m,n 
positive integers. It is easy to show that, in the absence of gap in the ground-state local-field distribution, the 
terms with i = j yield the dominant low-temperature corrections and vanish linearly with T. Terms with i =/= j 
yield higher-order 0(T 2 ) contributions. Suppose P({h*}) stands for the ground-state local-field probability 
distribution. For the terms x^x™ (i ^ j), we have 

/oo 
exp{-2mPh*a*)exp(-2n/3h*a*)P{h* 1 ,..,h* v )dh* 1 ..dh* v . (32) 
-oo 

The field variables h* k (k ^ i,j) may be integrated out, yielding the following expression 



x ,n x n = exp(-2m/3/i*Oexp(-2n^a-*)Pi J (/i* ) ^)d/i*dft* , (33) 
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/oo 
P(h$,..,h* v ) J] d K • ( 34 ) 

If the local field distribution P({h*}) has finite weight at the point hi = OVi, then the same holds for the 
two-sites probability Py(0, 0) so that we may expand this term around (0,0) in (33) obtaining thereby, 



x?x n = exp(-2m/3^V*)exp(-2n/3^V*)(P y (0,0) + ( -fiff) K + ( ) ^ + 0(ftJ/iJ)) dfcJd/$ 

(35) 

where 0(h*h*) denote higher-order terms at least quadratic in the fields. A simple saddle-point calculation (in 
the (3 — > oo limit) gives then, 

T 2 „ 

x?x? = — Py(0,0) + 0(T 3 ) . (36) 

The dominant terms in the limit T — > correspond to terms of the type a;™, which give 

— TPAO) , , 

71 

where P(0) is the value of the single-site probability distribution on the site i evaluated at h = 0, 

/oo 
i>(Jn,.. s /iv)$(/i i -fc*)IJdfc* . (38) 
-°° k^i 

This probability is not independent of the spin i, as our Hamiltonian can contain terms which introduce asym- 
metry between different sites. This is an important difference with respect to the computation of the previous 
section where the local field distribution ( |l3| ) was site independent. Actually, this independency was necessary 
in the "instructive example" to fully carry out the analytic computation of G and G c . The key point is that, 
at the level of one-spin excitations, low-temperature corrections to overlap averages are linear in T and Pi(0). 
According to expressions |3(]), ( |3l"| ) all sites are equivalent (inequivalence of sites enters only through the value 
of p(0)), so the only invariant term linear in P is ^ A(0). The numerator in (||) yields, 



-2 



16T^P(0), 



<9 2 ) 2 - <9 2 f = ^^r"' {V If + 0(T 2 ) . (39) 

To compute the overlap (q 4 ) we use the expression, 

W) = ^l(3V 2 ~2V+(6V-8)Y,Tv+ E T «") > ( 4 °) 

where Tijki = (aiaja^ai) 2 . Similarly as for the two-point correlation function ( p9| ) we obtain 

((Ti(Tjcrk<7l) = cr * cr j (T i o 'l0- ~ 2 ( x i + x i + x k + x i)) k, I all different) . (41) 

With the same assumptions as for the two-points function we obtain for the denominator in (|4|) 

W)-W) 2 = 16T y* (Q) (V - l) 2 + 0(T 2 ) , (42) 



which finally yields, 



G = I + 0(T) . (43) 



A similar calculation for G c yields G c — || + 0(T). 
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• Two-spin excitations 

Let us consider now excitations which involve only two different spins in the lattice (V(V — 1)/2 different type of 
excitations). In this calculation one-spin excitations are not included. It is easy to check that these excitations 
yield 0(T 2 ) corrections to the two-spin and four-spin correlations. Under the same conditions as before these 
are given by, 

(o-io-j) = a*a*(l - 4(x, + Xj ) x l) (* + j) ( 44 ) 
(ai<Tj(r k ai) = a*a*a* k a\{\~ ifa+Xj + x k + x{) ^ x m) (i, j, k, 1 all different) (45) 

m^i,j,k,l 

A saddle point calculation shows that corrections to the ground-state correlation functions are quadratic in 
T. Finite T corrections now depend on both Tlx] and XiXj for i ^ j. Now, for the quantity G we expect a 
dependence of both numerator and denominator on terms of the type Pi(Q)Pj(Q) as well as Pij(0, 0). They can 
enter in different forms, for instance J2i^j Aj(0, 0), Ei-Pi(O)) 2 or (J2i ^(0) 2 )- A universal value for G is not 
guaranteed anymore. In particular, supposing uncorrelated local fields (which in principle may not be true) and 
independency of the one-site probability distribution Pi (0) on the site i we obtain, after a simple but lengthy 
calculation, 

F) 2 - F) 2 = q^ (0)2 <Y 2) 2 (U - 1) + 0(T 3 ) (46) 
_ _ _ 2 = 64T^P(0) 2 (F _ 2)2{y + 0(r3) (4?) 

and their ratio yields G = | + 0(T) which is different than before. We stress again that the result 2/9 is not 
universal and will certainly not hold in the most general case. This calculation has been shown to stress how 
the 1/3 value is a fingerprint of the dominancy in the limit T — > of the one-spin excitations. 

• Higher- order excitations 

The generalisation to the most general case of i^-spin excitations is straightforward. Including only this class of 
excitations we obtain 0{T K ) corrections to correlations which involve any finite number of spins. This can be 
easily seen from the fact that any possible excitation of this type will involve the reversal of K different spins, 
each spin i contributing by a factor x\ = exp(— 2f3h*) to the correction. The simultaneous effect of all spins 
yields a product type Jlj=i x % w bich immediately gives (in the limit (3 — ■+ oo) the T K term. The numerator 
and denominator in G are of order T K with 0(T K+1 ) corrections. The final result for G for any value of K is 
not easy to compute and, as previously discussed, will depend on a larger number of invariants which involve 
different combinations of the terms Pi(Q), P^fi, 0), .., Pi lt ..i K (0, 0). 

When all possible excitations are treated together the calculation proceeds as before. The dominant contribution 
for OPF will always come from samples whose lowest excitations are one-spin excitations. Consequently, in the 
zero-temperature limit (for V finite) one-spin excitations dominate the correction to correlation functions proving 
our conjecture. Note that the result we are stating here is quite natural. OPF at very low temperatures are always 
dominated by those rare samples characterised by local fields f3h « 1 where one spin-excitations yield the largest 
contribution. From a numerical point of view this implies that more samples are needed to compute with a reasonable 
precision the values of G and G c as T goes down. This is because for T — > the effect from rare samples on OPF 
becomes more and more important. Let us stress again that the present derivation assumed that TV C 1, In the 
opposite limit or in an intermediate regime the result obviously does not hold. In that case, it may well happen that 
dominant contributions in OPF involve the reversal of a large number of spins (domain excitations) which, in the 
limit TV 1, may also involve the whole systemEil. 

The hypothesis of a unique-ground state is apparently in contradiction with the case in which there is time-reversal 
symmetry. Indeed all spin-correlations computed in this section are invariant under time-reversal symmetry and the 
present conclusions remain unchanged. The situation is certainly different in disordered systems with non-trivially 
degenerate ground states (for instance, finite-dimensional spin glasses with discrete couplings) where we expect that 
G(V,T) vanishes exponentially with 1/T like in the instructive example of the previous section. Again, in the other 
limit (finite tepaperature and V — > oo) the behavior of these degenerate models may completely change and G could 
be finite agairJla. 
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V. THE ID ISING SPIN GLASS 



In this section we present an analysis of the lD-Ising spin-glass model with free boundary conditions. We consider 
the folllowing hamiltonian: 



n 



N-l 

E 

i=l 



(48) 



where the couplings are randomly distributed according to the probability distribution P(J). Our aim is to obtain 
an analytic expression for G and A eqs. (|J) and (jit). As this model has the transition at T — 0, we expect that in the 
large volume limit G will go to zero except at T = 0, where G — 1/3. Moreover, we show that at zero temperature 
G = 1/3 for any finite system, although here the two limits (V — > oo and T — > 0) do not commute. In order to obtain 
an expression for the moment of the order parameter q we have computed the following object: 



( e yi) m , (49) 

where m is a positive integer and q is the overlap between two different configurations of spins, which is the generator 
of the moments of the overlap {q p Y ■ Once obtained this expression, by partial derivation respect to y we will obtain 
expressions for the expectation values of all the moments of q, such as: 



■7— r d n ( e yi). 



(50) 



In our computation we are only interested on the quantities: (q^Aq 4 ) and (q 2 ) 2 . Consequently we only need to 
compute d49| ) for m = 1,2. The former can be easily computed by (p0[). By doing some more work we can obtain an 
expression for (q 2 ) 2 : 



{q 2 ) 2 



d 4 {evi) d 4 {evi) 



dy 4 



dy 4 



(51) 



y=o 



where we have used the fact that in this model (q) = 0. 



A. The transfer matrix method 

For general m, ( ^9| ) can be computed through the transference matrix methodS We have to compute: 

n Y 2 (52) 

a=l 

where Z— 2 TJ i 2 cosh(/3 Jj) is the partition function of a ID chain, a is the index for each pair of replicas and we have 
m systems of two replicas. 

In order to perform the average over the disorder, we are interested in considering the transfer matrix associated 
to each point i, so that it contains all the dependence of the Jj. For a single pair of replicas this matrix reads: 

exp h/ CT ' r ' + ^ lT ' +1 + /3Ji( o-jo-j+i + TiT i+ i)\ 

V i ^V(a i ,n;* w ,n +1 ) = (acosh^)) 2 ■ (53) 

For general m our matrix associated to each point consists of the tensorial product of m matrices Vi- At this stage 
we are ready to perform the average over the disorder and for any i we have: 



T = Ti = Q9 Vi (54) 
i 
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Then our calculation is reduced to: 



(eviy 



1 , — , / .„ "i T l N-l / "n t n 

ly^e^^nv T e v^^2N 



so we must compute the trace of the product 



T N1 B 



(55) 



(56) 



where A is a 4m x 4m matrix, which is the tensorial product of m matrices, which contain the terms of the two edges 
which had fallen out in the symmetrization process, 



D 



l v zlzimil 
p y 2N 

2 2 



(57) 



The rest of the calculation is straightforward: In first place we have to diagonalise T, and obtain the set of eigenvalues 
and eigenvectors, so that in this new base we have 



N-l 



N-l 



, N-l 



\ 



\ 



(58) 



where the subindex A stands for the diagonalised matrix. We then have to obtain the change of base matrix M which 
expresses the new set of eigenvectors {A 1 } in terms of the old base {<J a }. We finally obtain: 



N-l 



{eyi) m = Tr MT\ M B. 



(59) 



We have to point out that the case m = 1 is easy to solve . However, the case m = 2 turns up to be more difficult 
as the diagonalization of V is not trivial. To compute (q 2 ) 2 one can always use the traditional method by using the 
fact that 



J^J tanh (3J P i ^ j 



(60) 



B. Results 

Here we report on the obtained results in the low-temperature limit and in the infinite- volume limit. The relevant 
quantities (q 2 ), (q 4 ) and (<7 2 ) 2 only depend on N, R and R 2 which have been introduced in section III, and whose low- 
temperature behavior is given by (^l|). At finite temperature, where R and R 2 <C 1 we obtain for the numerator 



and denominator in (Q): 
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4(1 + R)(R -R 2 ) ml 1 . 

numerator — = = ^ Oi—rr) , (61) 

iV 3 (i?-l) 3 (i? 2 -l) N 2 ' ' 

4(1 + R)(R 2 -R 2 ) , nl 1 . 

denominator = — — A '- + 0(—z) , (62) 

N 2 (R- 1) 3 {R 2 - 1) N 3 ' 

JV N 

where we have kept the lowest orders in 1/N and we have made the following aproximations linijv^oo R , R 2 — > 0. 
We see that in this limit G goes to zero as 1/N. However, if we take the low temperature limit (|l]), where A, A 2 w 1 
then we get the expressions 

AD(N 4 - l)T 

numerator = v 3 ' + 0{T 2 ) , (63) 
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denominator 



4D(N 4 



1)T 



157V 3 



C(T 2 ) 



(64) 



where D is given by D = 2P(0). This yields G = ± + 0(T), independently of the size of the system. A detailed 



computation up to second order in T gives us that in the large volume limit: 
fact for the parameter A, we get in the limit T — » 0: 



A = 



4L>(iV 4 - \)T 
45iV 3 



0{T 2 ) 



G = 4 — BTN, B being a constant. In 



(65) 



In figures (||) and @ we show G and A as a function of the temperature for a Gaussian distribution odf couplings 
P(J) = (2II) _1 / 2 ea;p(— J 2 /2). Note that the low-temperature correction to G and ( |S5| ) scale as IW when TV — > cxd 
reflecting the fact that as we get close to the transition point T = 0, the correlation length diverges as 1/T. We 
recover the desired result at T = 0, however we have to stress out that in this model both limits T — > and iV — » oo 
do not commute. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

T 

FIG. 5. Parameter G for the ID Ising spin glass for lengths N=2,4,8,16,32 (from above to below) 
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VI. THE SPHERICAL SHERRINGTON-KIRKPATRICK SPIN GLASS 



In this section we present some numerical simulations for the values of G and A in the Sherrington-Kirkpatrick (SK) 
spherical spin glass. This case is quite interesting because its low-temperature behavior corresponds to the second 
possibility mentioned in section II where OPF vanish (in the V — > oo limit) much slower than the paramagnetic 
example studied in the previous section. Correspondingly the study of OPF in this model turns out to be very 
complicated because the equilibrium solution is marginally stable. The model is defined bycJ 

H = — } j Jij<Ji<Tj , (66) 
where — oo < Oi < oo and the values of Oi satisfy the spherical global constraint Jw=i of = N. The. couplings have 



average zero and variance 1/N. The statics of this model can be solved with and without replicasEil. In the former 
case one finds a transition at a temperature T c = 1 where the Edwards- Anderson parameter is different from zero 
and equal to 1 — T. In the latter case the transition corresponds to a macroscopic condensation of spin configurations 
onto the eigenvector corresponding to the largest eigenvalue. In the replica framework it has been shownO that the 
replica symmetric solution is the only possible one within the Parisi scheme. Since OPF vanish, the computation of G 
requires the knowledge of finite-size corrections in the numerator and denominator in (j^) . A simple calculation reveals 
that the replica symmetric solution is marginally stable (the replicon eigenvalue vanishes everywhere below T c ) so the 
spin-glass susceptibility diverges. The situation is similar to what happens in the usual Sherrington-Kirkpatrick model 
with Ising spins. There the spin-glass susceptibility diverges proportionally to the volume while now the divergence 
is much slower (like N 1 / 3 ). This is so because in the present model OPF vanish like N~ 2 / 3 while in the original SK 
model OPF are finite. 

Again, to compute G we need to know the precise value of the amplitudes entering in the finite-size corrections in 
the parameters (q 2 ), (q 2 ) 2 , (q 4 )- It is well known that analytic calculations of finite-size corrections in spin glasses are 
extremely difficult, specially for the amplitudes which ape. the quantities we are interested in. For the SK model these 
amplitudes are partially known only for some quantities!!!. For the present case we will use theoretical considerations 
and numerical simulations to estimate the asymptotic behavior of the different overlaps. 

We have simulated model ( [66] ) with a Monte Carlo dynamics where a change of a randomly chosen spin is proposed 
ft ~~ * <J i + $ r i where 5 is a constant number typically of order 1 and is a random number uniformly distributed 
between —1/2 and 1/2. The value of 5 is chosen to have a reasonable acceptance rate. The value of all other spins is 
recalculated in order to satisfy the global spherical constraint. Moves are accepted according to the Glauber algorithm. 
Note that although we need to recalculate the value of all spins (changing them by multiplying by a normalization 
constant) the change in the energy can be simply calculated in a finite number of operations independent of N and 
simulations are as fast as with Ising spins. Our investigation has focused on small sizes, which reveal how G is a 
powerful tool to investigate phase transitions. The number of samples simulated are typically several thousands for 
very small sizes (N = 4, 6, 8, 12, 16) and several hundreds for larger ones (N — 24, 32, 40, 48, 64). Overlaps have been 
computed by collecting statistics over a large time window (tipically of order 10 5 Monte Carlo steps for each sample). 
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We have evaluated, {q 2 ) , {q 2 ) 2 , {q 4 ) for different sizes and temperatures. 

Figure shows the results for G. Note that already for the smallest sizes there is a crossing of the different curves. 
The crossing appears for values of T well above T c = 1 for the smallest sizes and moves to lower temperatures as 
the size increases converging to the expected value T c = 1. It is quite surprising that already for very small sizes the 
transition can be clearly seen. The crossing moves down in temperature as the sizes increase and already for several 
tens of spins converges to the correct value T — 1. As a comparison we show in figure ^ the behavior of the usual 
Binder parameter defined as 

*4(3-5) (67) 



2 



(Q 2 ) 



In this case the crossing point appears at low temperatures for small sizes and moves up very slowly as the size 
increases. But already for the largest sizes the crossing is still at T ~ 0.8 quite far f**Bi T = 1. A similar effect has 
been observed in simulations of the Sherrington-Kirkpatrick model with Ising sping!3'Q. These results indicate that 
a numerical study of thej-narameter G can be extremely useful to locate phase transitions in disordered systems by 



studying very small sizes 
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FIG. 7. Parameter G for the SK spherical spin glass (N = 4, 6, 8, 16, 24, 32, 48, 64 from below to above at low temperatures). 
The larges error bars are shown for the largest size N = 64. 
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FIG. 8. Binder cumulant B for the SK spherical spin glass (N = 4, 6, 8, 16, 24, 32, 48, 64 from below to above at low temper- 
atures). Error bars are now negligible. 



To analyze better the behavior of G at l ow temper atures we have tried to extrapolate G to the large N limit. Below 
T c we expect for all three quantities (q 2 ), (q 2 ) 2 , (q 4 ) the following finite-size corrections, 



(q 2 )\ (q 2 ) 2 , (q 4 ) = q% A + + - + + , 

with qEA = 1 — T. From these expressions we expect for G the following behavior, 

G = G X + J ^ + J ^ + 0(1/N) . 



(68) 



(69) 
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We have fitted the values of G to this expression with Goo, A, B as fitting parameters. The results and the fits 
are shown in figure [| The extrapolated values for the lowest temperatures T = 0.6, 0.7 are Goa{T = 0.6) = 0.34 ± 
0.2, (A(T = 0.6) = -0.71 ± 0.1 and B(T = 0.6) = 0.49 ± 0.13), G^T = 0.7) = 0.29 ± 0.2, (A(T = 0.7) = -0.66 ± 0.1 
and B(T = 0.7) = 0..49 ± 0.12). Within errors these are compatible with the value 1/3. Trying to have an estimate 
of Goo at higher temperatures is very difficult because critical effects arc strong. 
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FIG. 9. Fit function (p9|) to the G parameter for different sizes at T — 0.6 (above) and T = 07 (below). Extrapolations to 
N — > oo are compatible with G(V — * oo) = 1/3 in the spin-glass phase 

We must conclude that for this model the universal value 1/3 is well compatible with the data suggesting that 
this may be a generic result for a spin-glass phase. Still we should do more extensive simulations to reach a final 
conclusion. Although going to larger sizes at the lowest temperatures may be factible this will require much longer 
computational time. 

VII. OUTLOOK AND DISCUSSION 

In this paper we have investigated order parameter fluctuations (OPF) in spin glasses. In particular we have 
considered four different parameters: G, A for disconnected thermal averages and G c , A c for connected thermal aver- 
ages. Itrbas been recently shown that these r&edfils can be very useful to investigate phase transitions in disordered 
systems™ and several recent numerical worksQcTEl indeed support this conclusion. In this work we have concentrated 
our attention to obtain general results and to apply them to certain solvable cases where these can be explictly 
checked. 

We have demonstrated that for models with a unique ground state and no gap in the ground-state local field 
distribution (for instance, all discrete models with continuous disordered couplings taken from a distribution without 
gap) G and G c take the respective universal values G = 1/3, G c = 13/31 at zero temperature for any finite volume. 
This is consequence of the dominancy of one-spin excitations in OPF. For infinite volume this result still holds only 
in the regime where the limit T — > is taken before the limit V — > oo and fast enough such that TV — > 0. This result 
has then been checked calculating OPF in an instructive example without many body interactions and for the case 
of the one-dimensional Ising spin glass where explicit computations can be done using the transfer matrix method. 
All these good properties suggest that both parameters G, G c are ideal candidates to investigate phase transitions in 
disordered systems much alike the Binder cumulant is for ordered systems. 

The extension of this result to the other limit where V — > oo is taken before T — > or, more generally, the limit 
V — > oo for T finite is far from trivial. In this last case, G(V, T) is not volume independent anymore. So the question 
is whether G(F,T) converges in the large V limit to the universal temperature independent value 1/3. At finite 
temperatures there are different possible scenarios for the value of G. In case OPF arc finite in the V oo limit 
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stochastic stability arguments and replica equivalence suggest that G should be 1/3 everywhere in the spin-glass 
phase. Replica equivalence is a very generic property which, to our knowledge, has not been emphasized before in 
the present context and implies that the free energy of a replicated disordered system must be proportional to the 
number of replicas. Note that at zero temperature replica equivalence cannot be used because the limits V — > oo and 
T — » may not commute in that case. Actually, as we proved in section IV only for models with a unique ground 
state and absence of gap in the fields distribution, G takes the universal value 1/3 but vanishes (exponentially fast 
with 1/T) in the presence of a finite gap in that distribution. 

The other interesting case is when OPF vanish. And here we can offer only more speculative arguments. A possible 
scenario is that which distinguishes two possibilities depending whether, in the infinite-volume limit, OPF vanish like 
1/V or slower like 1/V a with a < 1. If OPF vanish like 1/V, G may take the value typical of a paramagnetic phase 
(for instance the case of the one-dimensional spin-glass model) or a temperature dependent value (the instructive 
example of section III). For these two solvable cases the parameter G is quite different. In the one-dimensional Ising 
spin glass we find G = |<5t,o while in the instructive example we find G = G(T) with G a monotonous decreasing 
function of T with G(0) = 1/3. The reason for these two different behaviors in a disordered phase may be adscribed 
to the fact that, in the first case, there is a critical point at T = while in the second there is no critical point at 
all. So G is a good indicator for a phase transition. But this observation must be taken with caution because the 
parameter G c shows a different behavior for the instructive example G c — -gf^T.o similar to the behavior of G in the 
one-dimensional spin glass. We expect the interesting behavior to be present in models where OPF vanish like 1/V" 
with a < 1. This class of models includes disordered systems where the replica symmetric solution is marginally 
stable and eventually finite-dimensional spin glasses if replica symmetry is not broken, a question still unsolvedc3. 
This case is much more subtle because replica equivalence cannot be used (nor probably the stochastic stability 
property) and finite-size corrections must be known. To address this question we have done a numerical study of the 
spherical Sherrington-Kirkpatrick spin glass. Two are the main outcomes: 1) The parameter G is an excellent tool 
to locate the spin-glass transition already for very small sizes (more precise than the usual Binder parameter) and 2) 
An infinite-volume numerical extrapolation (compatible with the expected form for the finite-size corrections) of the 
value of G in the spin-glass phase is well compatible with the value 1/3. 

shares the same p 

Before concluding we want to stress that, apart from their applicability to the study of spin-glass transitions, OPF 
are interesting quantities which deserve further investigation. The outcome of the proof in section IV is that OPF are 
much sensitive and rely completely on the effect of rare samples. Actually, rare samples are those which induce the 
largest OPF and fix the value of G to 1/3. A comprehensive study of rare events in disordered systems is still missing. 
Averaging of extensive quantities such as the replicated free energy in standard renormalization group approaches may 
wipe out a large number of effects such as those discussed here. Certainly more detailed investigations are needed to 
clarify the situation. Although a final theorem which resolves this problem may be at hand we think that the search 
for non-trivial counterexamples of the different possibilities discussed in this paper could be very useful. 
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